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SOME RESEARCH RELATED TO THE STRUCTURAL DYNAMICS OF LAUNCH VEHICLES 




By Harry L. Runyan 


NASA Langley Research Center 


SUMMARY 


\ 



This paper treats recent advances in the use of dynamic models for struc- 
tural dynamics and aeroelastic structures . Specifically, model tests of the 
Saturn I are discussed as ¥611 as tests of the l/lO- and l/40- scaled models of 
the Saturn V. The problem of the effect of ground winds on launch vehicles is 
discussed, including a description of the model program and the results of some 
tests. 


Also, the problem of the high-frequency buffet response of a launch vehi- 
cle adapter is treated, including a comparison of flight measurement with 
calculations. 


INTRODUCTION 


£ 



The purpose of this presentation is to concentrate on structural dynamics 
and loads which are of importance in the design and operation of launch vehi- 
cles from the time of lift-off to the time at which the vehicle leaves the sen- 
sible atmosphere. This time interval is relatively short, being of the order 
of just several minutes. 

What are some of the problems which are of concern to the dynamicist in 
designing and operating a launch vehicle? A listing of these problem areas is 
given in figure 1. On top of the figure is shown "input" loading, which is 
imposed on a "system" - the launch vehicle, and the "output" load, in terms of 
stresses and acceleration. 

Buffet and acoustic loads are closely related as can be seen by the gen- 
eral similarity of the outputs. Winds, both ground and aloft, can induce 
bending moments which can involve such severe loads that certain portions of 
the vehicle must be designed to accommodate these loads. 

Fuel slosh in general is a control problem; however, loads can be induced 
into the vehicle if the liquid oscillations are not kept to a minimum. Control 
loads during programed maneuvers are sometimes severe and, in particular, 
excessive loads from maneuvers should not be applied during periods of high 
loading from other sources. Recently, response and loads in the longitudinal 
direction have become of increasing concern. Since 90 percent of the weight of 
a typical liquid-propellant vehicle may be liquid, the coupling of the liquid 
system with the structure and engine is possible. During lift-off, when the 
hold-down clamps are released, a large dynamic response is found and is the 
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principal design consideration for the bottom of the tanks. Another dynamic 
phenomenon which has recently been found on two large vehicles is the coupling" 
of the engine, turbopump, structure, and fluid components into a system that , 
can be unstable. This has been referred to as the POGO problem, since the vehi- 
cle responds in a longitudinal direction, in the manner of a boy who is 
bouncing on a pogo stick. 

One of the unfortunate aspects of this loading picture is the simultaneity 
of the occurrence of some of these loads. In figure 2, a plot of various 
loading sources against time of flights is shown, where the black area indi- 
cates the probable time at which the particular event will occur. The point 
here is that it appears rather black between the transonic region and the maxi- 
mum dynamic pressure region for most of the loading conditions. These combined 
loads present a rather complex interaction problem which, at the present time, 
is not entirely understood. 

To illustrate this loading picture in more detail, on figure 3 is plotted 
bending moment versus vehicle station. The curve labeled "capability" repre- 
sents the bending-moment capability remaining after such known loading factors 
as internal pressure, acceleration, and drag have been accounted for. Thus, 
the loads from the remaining sources can be summed to ensure that their total 
does not exceed the capability curve. Of particular note is the large incre- 
ment in bending moment due to winds aloft. At certain times the winds may be 
very severe, and may be as high as 100 meters per second or, in rare instances, 
even higher. 

To illustrate this wind picture lines of constant wind velocities for the 
Northern Hemisphere are presented (isotacks) as obtained from reference 1. A 
typical result is shown in figure This figure represents the average wind 
over a 3-month period for an altitude of 11.8 km. Statistically, these veloc- 
ities will be exceeded only 1 percent of the time. Note the low wind velocity 
found both on the Equator and in the Arctic region. Of particular Interest, 
however, since this meeting was held, in Tokyo under the auspices of the Japanese 
Rocket Society, is the high wind region over Japan. The 99 percent wind veloc- 
ity is about 90 meters per second, which is the highest wind speed of any area 
in the Northern Hemisphere. 


DYNAMIC MODELING 


A basic ingredient in the solution of most structural dynamic problems is 
a knowledge of the vibration modes. Eor complex cases, theory is not adequate 
and the use of dynamic models is mandatory. The next section of the paper will 
be concerned with specific examples of dynamic model testing, concentrating on 
the problem of the gross vehicle vibratory motion. 
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l/5-Scale Vibration Model of Saturn 

. Description , - The l/5-scale Saturn model described in reference 2 was 
constructed for tbe study of lateral tending vibrations j therefore, the impor- 
tant parameters to be scaled were the mass-stiffness ratios. Such parts as 
aerodynamic fairing and fuel piping were not scaled since they did not contrib- 
ute to stiffness j however, lead ballast weights were used to simulate their 
mass. Furthermore, the vibrations of principal interest were the overall vehi- 
cle modes, so local pauex were not scaled. Other factors, such as 

fuel sloshing, were considered to have a secondary effect on the overall vehi- 
cle vibration and so were not scaled. The model is shown on the left in fig- 
ure 5j and the full-scale vehicle is shown on the right in a vibration test 
tower . 

The type of scaling chosen for the Saturn model was a component -by- 
component uniform reduction of dimensions to l/5 of the full-scale values, 
using the same materials as the full scale. This replica type of scaling was 
chosen because of the structural complexity of the Saturn booster with the 
resulting difficulty of determining accurate equivalent stiffness and mass 
properties for the multiple-beam truss-work assemblies incorporated in the 
vehicle. 

Results . - Some results of the model vibration test are shown in figures 6 
and 7* In figure 6, the first or fundamental mode shapes and frequencies 
of the model and full-scale vehicles are shown with the vehicle ballasted to 
simulate the maximum dynamic-pressure weight condition. This figure shows 
almost exact agreement between model resonant frequency, when adjusted by the 
scale factor, and the full-scale resonant frequency. The mode shapes of the 
full -scale vehicle and the model are in good agreement as can be seen from com- 
parison of the circles with the square symbols. The booster outer tanks, which 
are free to deflect independently of the center tank except at their ends, are 
shown in this figure to have the same deflection as the center tank, and this 
mode has the appearance of the more conventional bending modes obtained from 
nonclustered vehicles. 

In contrast, the second vibration mode, shown in figure 7> shows one of 
the unusual vibration modes associated with the clustered arrangement of the 
Saturn booster. There is about a 10-percent difference in frequency between 
model and full scale in this mode, and comparison of the circles with the 
squares shows fairly good agreement of model with full-scale mode shape. The 
typical outer tank indicated by the flagged symbols is seen to deflect in the 
opposite direction as the center tank. Cross section A-A at the midsection of 
the booster shows that when the center tank deflects upward as indicated by 
the arrow, the outer tanks are deflecting independently in the downward direc- 
tion. This unusual mode shape where the outer tanks deflect in the opposite 
direction from the center tank is associated with the clustered arrangement of 
the booster tanks and has been termed a "cluster mode." 
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Saturn V Dynamic Models 1 

Figure 8 illustrates the present modeling project, which is a l/lO-scale 
replica model of the Saturn V. The model is 36 feet high, including the Apoll'o 
payload. The model will duplicate all of these essential details up to the 
Apollo payload. The general type of construction of the full scale including 
skin stringer, waffle pattern, integrally milled stringer, and corrugated skin 
is being followed. As an example of the construction detail adhered to in the 
model, figure 9 illustrates the interior of the first-stage tankage. Note in 
particular the detail in the fuel slosh baffles which are the thin corrugated 
rings located on the circumference of the tank. The three corrugated pipes are 
three of the five tunnels through which the fuel lines pass through the bottom 
tank to the five individual pumps for each engine. 

In addition to the l/10-scale model, a l/40-scale model is being construc- 
ted, which is shown in figure 10 and which is about 9 feet high. The modeling 
concept for this case involved only the duplication of the mass and stiffnesses, 
and does not follow the replica concept. 

A sampling of the results of vibration testing of the two scaled models 
along with theoretical predictions are given in figure 11 taken from refer- 
ence 3* Here the lateral vibration frequencies for the first three modes are 
plotted as a function of the first-stage propellant load, with the upper stages 
remaining full. The experimental results are plotted as open symbols and the 
analytical results are shown as solid or dashed lines. Examination of the 
experimental data indicates excellent agreement between the two models, in 
spite of the fact that one model was essentially an exact replica of the full 
scale (l/lO), whereas the other (l/4o) was designed solely on the basis of an 
equivalent stiffness and mean distribution. 

The analytical results are interesting from the standpoint of the pretest 
results. The pretest results were based on estimation of the stiffness and 
mass from early design drawings. The post-test results were obtained after the 
model had been constructed, parts weighted, and better and more detailed stiff- 
ness estimates made. The increase in accuracy of the analytical results is 
apparent, and may be considered acceptable for design purposes. 


BUFFETING 


The term "buffeting" is really not very well defined and precise. For the 
purposes of this paper, the term "buffet" will refer to aerodynamic flows 
involving departure from potential flow with consequent flow breakdown and 
separation which can result in structural response of the vehicle. Obviously, 
this is a rather broad definition and is symptomatic of the wide spectra of 
buffet problems. 

The occurrence of buffeting on launch vehicles has been a surprise to 
many people; the long cylindrical shape of a typical launch vehicle did not 
appear at first to offer any substantial problem. However, certain vehicle 
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failures in the early days of the space program were diagnosed to have been 
‘due to possible flow breakdown, resulting in structural failure. (See ref. 4.) 


As an illustration of the response of a structure to high-frequency buffet, 
scane work is in progress on attempts to calculate the random response of a com- 
plex structure such as a launch -vehicle adapter for comparison with results 
measured in flight. For this example, the comparison of the calculated and 
measured response of the adapter used on an early Mercury flight will be shown. 
The adapter is a ring and a Lriugcr stiffened shell structure which connected 
the spacecraft to the launch vehicle. The Mercury had an abort escape tower 
containing a solid-propellant rocket projecting ahead of the spacecraft. The 
combination of the rocket, the open tower, and the shape of the Mercury capsule 
resulted in some rather severe buffeting loads on the adapter immediately below 
the capsule. Measurements of the acceleration response of a point on the 
adapter were made in flight, and efforts have since been made to estimate the 
response spectra. The following results were obtained under contract to NASA- 
Langley Research Center by the McDonnell Aircraft Corporation, and are reported 
in reference 5* 


The fundamental equation for the cross-power spectral density is 


= |0alJ[sq J (^|j^ 

where 

Sq J (ia>) = f 2 f f 2 f [^(xjejX'je'jioJfSj^eJfSjU' ^'MxMx' )dfi < 1 x 6 ' dx"! (l) 
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Sp is the cross-power spectrum of the pressures defined by the cross- 
correlated function 


Rp(x,0,x* ,9' ,t) 


lim 
T -» oo 


1 _ 

2T 



p(x,0,t + r)p(x’ ,0 * ,t)4t 


where 


Sp(x,9,x' ,0’ ,±to) 



e _iujr Rp(x,0 ^x'^O'jt) cLt 


( 2 ) 


( 3 ) 


The basic and difficult problem is the accurate and practical integration 
of equation (l). There are two difficult quantities to determine. First, is 
the aerodynamic input characterized by the cross-spectra Sp. No analytical 
methods are available for its prediction in complex cases such as the Mercury 
configuration. For boundary -layer phenomena, estimates can be made based on 
many past experiments. For this case, however, resort must be made to measured 
quantities. The input pressure spectra were estimated from wind-tunnel meas- 
urements utilizing a 7 -percent scaled rigid model of the Mercury-Atlas config- 
uration including the escape tower. The vibration modes of the adapter were 
determined analytically, but the frequencies were adjusted to agree with 
results of ground vibration tests of a full-scale adapter. One of the signifi- 
cant factors in such an analysis is a good estimate of the structural damping, 
since the response is a direct function of the reciprocal of the structural 
damping at a resonant frequency. Damping is a very difficult and tenuous 
quantity to determine; many factors influence its value such as temperature, 
quality of joints, end condition, amplitude, etc. For this case, a nominal 
value of 2-percent critical damping was assumed. The integration required in 
equation (6) was accomplished by dividing the structure into 112 subareas, 
adjusting the properties across each subarea and proceeding with the necessary 
integrations as described in detail in reference 5* 

The comparisons of flight and calculation are shown in the next two fig- 
ures. In figure 12 is plotted the acceleration in g*s rms (0 to 600 cps) 
versus Mach number of a point on the adapter. The curve represents the flight 
data; the points represent calculated data. The agreement is quite good con- 
sidering the numerous approximations made in the analysis. Actually, what is 
being compared here is the total area under each acceleration spectrum at a 
given Mach number, since the area under a spectrum represents the mean square 
of the quantity. Information about the frequency content must be obtained 
from the spectrum itself. 

In the next figure (fig. 13) is shown the acceleration spectrum in g^/cps 
versus frequency in cycles per second for Mach number of unity. The dashed 
line is the result of the analysis, and the solid line is the value obtained 
from flight. Here the comparison is not as good as in the previous figure; 
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'however, the results are still considered excellent. For instance, the magni- 
• tude of the peaks is very similar, and when one considers such important fac- 
tors as the variability of structural damping, the effects of acceleration 
*Loads on the adapter as well as bending loads due to winds aloft, all of which 
were neglected in the analysis, the correlation is very good. Examination of 
the maximum response peaks indicates a shift in frequency. Four large peaks 
are apparent in each set of data with the analytical peaks occurring at a 
lower frequency. 

Thus, the picture for the estimation of the random response of a complex 
structure seems to be brightening, providing such factors as structural damping, 
vibration mode shapes, and the random aerodynamic input are adequately and 
carefully determined. 


RESPONSE TO GROUND WINDS 


For a final topic, another case of structural response to random disturb- 
ances will be discussed, namely, the response of a launch vehicle (as it stands 
on its launch pad) to ground winds. This is really an age-old problem; civil 
engineers have been concerned with the oscillation of smokestacks, buildings, 
and suspension bridges for a number of years. In spite of all this background 
work, it is still not possible to analytically predict the motion of a struc- 
ture to ground winds. For large launch vehicles, this is a problem which must 
be solved. These wind loads create problems in structural strength, guidance 
alinement, and clearance between adjacent structures. A more thorough discus- 
sion of this problem is given in reference 6. 

In order to obtain immediate answers for engineering purposes, resort to 
the testing of dynamic models in wind tunnels and the measurement of the 
response is necessary. But before discussion of the program, what is the basic 
problem? In figure l4 are shown schematically a launch vehicle, a tower, and 
a wind profile. The wind profile (on the right) can be considered as a steady 
wind having a profile varying in the vertical height upon which may be super- 
imposed turbulence. The wind may or may not pass near or through an adjacent 
launch tower before impinging on the vehicle, but the presence of the tower 
is vital in studying this problem. The steady -wind component may impose both 
static and dynamic loads on the vehicle. The static load response is primarily 
in the direction of the mean wind, and the dynamic response which is associated 
with vorticity shed from the vehicle is largest in the direction perpendicular 
to the wind. These dynamic loads, which may result from either a random forced 
response or a periodic self -excited response, may be several times larger than 
the static drag loads. 

To provide an idea of the extent of the effort on ground winds , figure 15 
shows nine different configurations which have been Investigated in the Langley 
Research Center l6-Foot Transonic Dynamics Tunnel, ranging from a 15-percent- 
scale Scout vehicle to a 3 - P ercen t-scaled Saturn V. Full-scale Reynolds number 
is s imula ted for all vehicles except the 3 -percent Saturn V which was tested at 
about one -third full-scale Reynolds number. 
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Some interesting results obtained on the 3-percent model of the Saturn V ” 
are shown in figure 1 6 . The maximum lateral oscillatory bending moment is 
plotted against the Strouhal number U/fD, where U is the mean steady veloc- 
ity , f is the fundamental bending frequency, and D is the diameter of the " 
base. 


The response is given for three values of damping . For the lowest damping 
case, a very large response was found over a very limited range of Strouhal 
number. As the damping was increased, this peak completely disappeared. It 
is believed that this peak response is a case of self -excited motion, as in 
flutter, whereas the response for the higher damping ratio is simply the 
response of a low damped system to random inputs. The typical response 
record is shown, the large peak response being of the nature of a pure sine 
wave at the fundamental frequency. The response for the two higher damping 
ratios appears as the response of a single -degree -of -freedom system to a random 
input. Research is continuing to better define this large response. 


CONCLUDING REMARKS 


The fundamental ingredient in the solution of structural dynamics problems 
is a knowledge of the vibration characteristics of the structure. A dynamic 
model program was described, which provides for input data for both direct 
design analysis as well as for correlation of theory and experiments. Two 
examples of problem areas were discussed in which the random nature of the 
input is an essential part of the problem. These were the buffet response of a 
localized area of a launch vehicle and the response of a vehicle to ground 
winds . 
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Figure 2.- Loading conditions versus time of flight. 



Figure 3 • - Bending-moment capability and applied loads. 






Figure 5«- Saturn vibration-test vehicles. 
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Figure 8.- l/lO-scale model of Saturn V. 





Figure 9 .- Internal view of Saturn V model fuel tank 





Figure 10.- l/^O-scale model of Saturn V. 
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Figure 12.- Comparison of rms acceleration. 
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Figure 13 . - Acceleration spectra. 
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Figure 14.- Load conditions caused by ground winds. 
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Figure 15*- Some dynamic models used in ground winds research. 





Figure 1 6 .- Effect of damping on Saturn V model response 



